Secondary instabilities of hexagons: 
a bifurcation analysis of experimentally 
observed Faraday wave patterns 
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Abstract. We examine three experimental observations of Faraday waves gen- 
erated by two-frequency forcing, in which a primary hexagonal pattern be- 
comes unstable to three different superlattice patterns. We use the symmetry- 
based approach developed by Tse et al. [1] to analyse the bifurcations involved 
in creating the three new patterns. Each of the three examples reveals a dif- 
ferent situation that can arise in the theoretical analysis. 



1. Introduction 

The classic Faraday wave experiment consists of a horizontal layer of fluid that 
spontaneously develops a pattern of standing waves on its surface as it is driven by 
vertical oscillation with amplitude exceeding a critical value. Recent experiments 
have revealed a wide variety of complex patterns, particularly in the large aspect 
ratio regime and with a forcing function containing two commensurate frequen- 
cies [2, 3, 4]. Transitions from the flat surface to a primary, spatially periodic, 
pattern can be studied using equivariant bifurcation theory [5] . These group theo- 
retic techniques may also be applied to secondary spatial period-multiplying tran- 
sitions to patterns with two distinct spatial scales (so called superlattice patterns) 
as demonstrated by Tse et al. [1]. 

We apply the method of Tse et al. [1] to the analysis of three superlattice pat- 
terns observed when secondary subharmonic instabilities destroy the basic hexag- 
onal standing wave pattern in two-frequency Faraday wave experiments. We can 
make use not only of the general symmetry-based approach from [1] but also of 
many of the detailed results. The reason for this is that in their paper, Tse et al. 
considered instabilities of hexagonal patterns that broke the translation symmetry 
of the hexagons, but that remained periodic in a larger hexagonal domain com- 
prising twelve of the original hexagons. The instabilities under consideration here 
satisfy exactly the same conditions (though in fact they remain periodic in smaller 
domains as well). 

We begin by specifying the coordinate system and symmetries we will use 
in section 2, then describe the symmetries of the three experimental patterns in 
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Figure 1. The coordinate system and certain elements of the 
symmetry group F. The origin of the coordinate system is at the 
centre of the diagram, and the point (1, 0) is indicated. The small 
hexagons represent the primary pattern, which is invariant under 
reflections {k^ and Ky), 60° rotations (p) and translations (ti 
and T2). The secondary patterns are all periodic in the larger 
hexagonal box. The three corner points labelled with solid circles 
are identified through the assumed periodicity. 



section 3. In section 4, we apply Tse et al.'s method of analysis to these three 
patterns, and present normal forms and stability calculations in section 5. We 
conclude in section 6. 



2. Coordinates and symmetries 

The primary pattern is made up of regular hexagons, which are invariant under 
the group Dq (made up of 60° rotations and reflections) combined with translation 
from one hexagon to the next (see figure 1). Tse et al. [1] studied experimental 
patterns reported in [6], which had the feature that after the secondary instability, 
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the pattern remained periodic in the larger hexagonal box in figure 1. The 144- 
element spatial symmetry group of the primary hexagonal pattern within this box 
is r, generated by the following reflection Kx, rotation p and translations ti and T2: 

( 3 \f2>\ 

Kx : {x,y) ^ {-x,y) n : {x,y) ^ {x,y) + { — ] (1) 



2' 2 



We also define Ky = k^P^, and note the following identities: 

= 4 = = ^1 = ^2 = = identity, (3) 

pKx = KxP^, TiKx = KxtIt2, T2Kx = KxT2, (4) 

TiP = PtIt2, T2P = pTl, nT2=T2Tl. (5) 

The time translation tt advances time by one period T of the forcing function, 

which is the same as the temporal period of the hexagonal pattern. This time 
translation is combined with the spatial symmetries above to give spatio-temporal 
symmetries. 



3. Experimental patterns 

The three experimentally observed patterns are shown in figure 2(a-c), visualised 
using the techniques described in [7]. Patterns (a) and (b) are both obtained 
using Dow-Corning silicone oil with viscosity 47cSt and layer depth 0.35 cm, while 
pattern (e) was foimd using a 23 eSt oil layer of depth 0.155 cm. All three patterns 
are obtained with forcing function containing two frequencies in the ratio 2 ; 3; 
pattern (a) is found with frequencies 50 and 75 Hz, pattern (b) with frequencies 
70 and 105 Hz, and pattern (c) with 40 and 60 Hz driving frequencies. Pattern 
(c) was reported previously in [7]. Typically, the secondary bifurcations occur at 
forcing amplitudes between 10 and 50% larger than the critical acceleration for 
the primary hexagonal state. Further experimental details can be found in [7, 8]. 

For the purposes of the analysis, we consider the idealised versions of these 
experimental patterns, shown in figure 2(d-f). The first pattern in figure 2(a,d) 
retains the Dq symmetry of the original hexagons but breaks certain translation 
symmetries. It is periodic in the medium-sized dashed hexagon in figure 2(d), 
which implies that the pattern is invariant under the translations rf and T1T2. 
It has no spatio-temporal symmetries. The second pattern is similar, although 
it possesses only triangular (D3) symmetry instantaneously. Moreover, it has the 
spatio-temporal symmetry given by a 60° rotation combined with advance in time 
by one period T of the forcing, as in figure 2(e,g). In fact, this spatio-temporal 
symmetry was first suggested by the analysis below, and found to be consistent 
with the experimental observations. The third pattern in figure 2(c,f) is quite 
different: the dark lozenges in figure 2(f) represent the enlarged gaps between 
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(g) (h) 




Figure 2. Experimental and idealised secondary patterns, (a- 
c) Experimental patterns, visualised from above, (d-f) Idealised 
versions of (a-c). (g-h) patterns (e-f) but seen one forcing period T 
later. The idealisations arc all rotated by about 30° compared with 
the experimental pictures. 
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Table 1. Character table of the group F, taken from Tse et al, 
with corrections. A representative element is shown on the second 
line for each conjugacy class (see also figure 3), and the number 
of elements in the class is on the third row. The next fifteen rows 
give the characters associated with each conjugacy class for each 
of the fifteen representations. 



the hexagons in figure 2(c). The pattern is periodic in the medium-sized dashed 
hexagon in figure 2(f), so is invariant under translations and = rf. It is 

also invariant under the group of symmetries of a rectangle D2, and possesses the 
spatio-temporal symmetry of the translation T2 combined with advance in time by 
one period T of the forcing, as in figure 2(f,h). 

Using the information above, we write down the instantaneous (spatial) sym- 
metry groups of the three patterns from figure 2(a-c) in terms of their generators: 

Sa = (Ka;,P,Ti ,Tir2), Efc (Kj; , , T^* , riT2) , T,^ ^ {Kx , KyT2 , T^) ■ (6) 

These groups are of order 48, 24 and 12 respectively. For the full spatio-temporal 
symmetry groups, we would also include prr in the generators of Sb, and T2Tt 
in the generators of Ec, but initially we will work with the spatial symmetry 
groups. The reason for this is that the instantaneous (spatial) symmetries can 
be determined reliably from a single experimental image, while extracting spatio- 
temporal symmetries from the experimental data is more involved. 

Each of the three instabilities that generates the three different patterns will 
be associated with a set of marginally stable eigenfunctions; the new pattern, at 
least near onset, can be thought of as (approximately) a linear combination of 
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these marginal eigenfunctions and the original hexagonal pattern. Which linear 
superpositions are consistent with the nonlinearity inherent in the pattern for- 
mation process is determined by our bifurcation analysis. The symmetries in F 
all leave the primary hexagonal pattern unchanged, so they must send marginal 
eigenfunctions onto linear combinations of marginal eigenfunctions, which induces 
an action on the amplitudes of these functions. In other words, if there are n mar- 
ginal eigenfunctions /i, . . . , /„, with n amplitudes a = (ai, . . . , a„) G M", each 
element 7 e F sends a to R-yS., where the set of n x n orthogonal matrices Ry forms 
a representation Rr of the group F. For subharmonic instabilities of the type of 
interest here, this will generically be an irreducible representation (irrep) [5] . Tse 
et al. [1] have computed all the irrcps of the group F; the character table of these 
representations is reproduced in table 1. Recall that the character of a group cl- 
ement 7 in a representation is the trace of the matrix R^, and that conjugate 
elements (which form a conjugacy class) have the same characters. 

Once the representation associated with each of the three transitions is iden- 
tified, we can write down the normal form, work out what other patterns are 
created in the same bifurcation, and compute stability of the patterns in terms of 
the normal form coefficients. 



4. Method 

The first task is to identify which representation is relevant for each bifurcation. 

Tsc et al. [1] outlined a two-stage method to accomplish this. First, any symmetry 
element that is represented by the identity matrix in a particular representation 
must appear in the symmetry group of every branch of solutions created in a bifur- 
cation with that representation. This can be used to eliminate from consideration 
any representation that has an element with character equal to the character of 
the identity that does not appear in the symmetry group of the observed pattern. 
Second, we make use of the trace formula from [5], which gives the dimension of 
the subspace of M" that is fixed by a particular isotropy subgroup S of F with 
representation given by the matrices Rr- 



where |S| is the number of elements in S. Specifically, we use the trace formula 

to eliminate those representations for which the spatial symmetry group of the 
pattern fixes a zero-dimensional subspace (implying that the subgroup is not an 
isotropy subgroup); only the remaining representations need be examined in more 
detail. 

We proceed by first counting the number of elements in each conjugacy class 
for each of the symmetry groups Sa, and Sc. Figure 3 shows representative 
elements from each class and is helpful for this categorization. The result of this 




(7) 




(g) K^Ti (12) (h) n^T2 (12) (i) K^rf (6) 




(j) KyTf (18) (k) p (24) (1) p2 (8) 




Figure 3. The 15 conjiigacy classes of T. One clement from and 
the number of elements in each class are indicated. The letters 
(a)-(o) correspond to the columns of table 1. 
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is: Tia contains: 

a : 1, 6 : 6, c : 6, / : 3, i : 6, j : 6, A; : 8, Z : 8, m : 1, o : 3 (8) 
(that is, one element from class a, six from class b etc.); S(, contains: 

a : 1, 6 : 6, / : 3, i : 6, / : 8; (9) 

and Eg contains: 

a : 1, 6 : 1, c : 3, e : 2, /i : 2, o : 3. (10) 

The clement does not appear in the symmetry groups of patterns (a) and (b), 
which eliminates representations 1-6 and 9-12 (since is represented by the 
identity matrix in all these: see table 1). Similarly, tiT2 in class / and do not 
appear in Ec, which eliminates representations 1-9, 11 and 13 from consideration 
for that bifurcation problem. 

Next, by applying (7), we find that pattern (a) has a non-zero dimensional 
fixed point subspace only in representation 7, as does pattern (b). The spatial 
symmetry group of pattern (a) fixes a one-dimensional subspace, and that of pat- 
tern (b) fixes a two-dimensional subspace. Pattern (c) has a one-dimensional fixed 
point subspace in representations 10 and 12, and zero in other representations. 

We are therefore faced with three different situations: the spatial symmetry 
group Sa fixes a one-dimensional subspace in representation 7, so we expect by 
the Equivariant Branching Lemma (see [5]) that such a pattern will generically be 
found in a bifurcation problem with that representation. 

Pattern (b), on the other hand, has a spatial symmetry group that fixes a two- 
dimensional subspace. However, we must take into account that the pattern arises 
in a subharmonic (period-doubling) instability, and extend the groups F and Sf, to 
the spatio-temporal symmetry groups that arise by including time translations. We 
may then show that the spatio-temporal symmetry group of pattern (b) fixes a one- 
dimensional subspace, and so also arises generically in a subharmonic bifurcation 
with representation 7. This is the same representation as with pattern (a), obtained 
for similar experimental parameter values. Extending to include the subharmonic 
nature of the instability does not affect the branching of pattern (a). 

The third situation arises with pattern (c), which on symmetry arguments 
alone could be associated with either representation 10 or representation 12. In- 
cluding information about the spatio-temporal symmetry of the pattern does not 
distinguish between these two representations. However, information on the Fourier 
transform of the pattern does allow a choice to be made between the two possibili- 
ties; in order to show this, we first need to work out which combinations of Fourier 
modes are associated with each pattern. 

It is useful to have sample Fourier modes for the basic hexagonal pattern: 

/o(ar,y) = cos2^(^y^ + cos27r (^-| + -^^ + cos27r (^-| - -^^ , (11) 

with wavevector of length as well as sample Fourier modes for representa- 
tions 7, 10 and 12. The method described by Tse et al. [1] yields Fourier functions 
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that would be included in the eigenfunctions associated with representation 7; 
representative functions with the shortest wavevectors include: 

/,(x.v) = C0S2. (I + JL,) +c„2. (I - ^) +™2, (^) (12) 

= ™2, (I + + ,to2. (-1 + + si„2, (-^) , (13) 

which is made up of wavevectors of length equal to of that of the basic hexagonal 
pattern. Eigenfunctions for representation 10 are made up of Fourier functions that 

include: 

with wavevector of length 5 the fundamental; and representation 12 has: 

with wavevector of length the fundamental. In each case, we have chosen the 
Fourier modes with the shortest wavevectors, as these are easiest to identify in an 
experimental Fourier transform. 

The images of the Fourier transform of pattern (c) in [7] show that the mode 
created in the instability contains wavevectors that are a factor of 2 shorter than 
the shortest in the basic hexagonal pattern, which is consistent with representa- 
tion 10 but not 12. In this way, information about the power spectrum of the 
pattern is necessary to supplement the arguments based entirely on symmetries 
and to distinguish between the two choices. 



5. Normal forms 



Using the functions specified above as a basis for representations 7 and 10, the 
matrices that generate the two relevant representations are, for representation 7: 



Rk^ — h, 



Rp — 







_1 
1 

2 



2 
_ 1 
2 



Rtq. — Rt, 



where /„ is the n x n identity matrix; and for representation 10: 



Rko: — 



Rti — 





1 


-1 









-1 

0" 



-1 



Rn 



Rr 





1 



-1 





-1 



1 



Rtt — ~^25 



(16) 



(17) 



Rr^ = -Is. (18) 
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The perturbation amplitude at time j + 1 times the forcing period, given the 
perturbation at time j, is given by aj+i = f (aj), where the equivariance condition 
amounts to Rjf{a) = f(iiya) for all 7 e F. Using this, we can determine the 
relevant normal form associated with these two representations: 

zj+i = -(1 + ^i)zj + P\zj\'^Zj + Q\zj\'^Zj + Rz^ (19) 

for representation 7 (truncated at quintic order), where the two amplitudes of /i 
and /2 in (12-13) are the real and imaginary parts of z, and P, Q and R are real 
constants. For representation 10 we truncate at cubic order and obtain: 

aj+i = - (1 + lJ.)aj + Pa^j + Q{a] + b] + c])aj , (20) 
6,+i = -(1 + ii)hj + Pb] + Q{a''j + 6| + (^j)bj, (21) 
Cj+i = -(1 + fj.)cj + Pc] + Q{a1 + h] + c])cj, (22) 

where P and Q are (different) real constants. In these two sets of equations, n repre- 
sents the bifurcation parameter. The —1 Floquet multipliers at /U = arise because 
these are subharmonic bifurcations. In representation 7, equivariance with respect 

to -Rxt ~ — 

I2 is a normal form symmetry, so even terms up to any order can be 
removed from (19) by coordinate transformations [9]. With representation 10, the 
matrix — J3 = i?^ appears as a spatial symmetry, so the normal form symmetry is 
in fact exact, and every solution branch has the spatio-temporal symmetry ttP"^, 
a rotation by 180° followed by time-translation by one period. 

The patterns arc neutrally stable with respect to translations in the two 
horizontal directions, and so also have two Floquet multipliers equal to 1 associated 
with translation modes. We have neglected these as all the patterns we find are 
pinned by reflection symmetries that prohibit drifting. 

The final stages are to determine the solutions that are created in each of 
these bifurcations, their symmetry and stability properties, and to compare these 
with experimental observations. 

The first normal form (19) generically has two types of period-two points, 
found by solving f{z) = —z: 

The first of these has exactly the symmetry group of pattern (a), with no spatio- 
temporal symmetries, while the second has exactly the spatial symmetry group Sf, 
of pattern (b), as well as spatio-temporal symmetries generated by prx- Recon- 
structions of these two are shown in figure 4(a) for pattern (a) and figure 4(b,c) for 
pattern (b), using the Fourier fimctions from above. Linearising the normal form 
about these two period-two points readily yields stability information: if P > 0, 
then both patterns are supercritical but only one is stable, while if P < 0, both 
are subcritical and neither is stable. 



Secondary instabilities of hexagonal Faraday waves 




Figure 4. Reconstructed patterns from the two solutions that 
arise in representation 7, using the Fourier functions (12-13) 
added to a function of the form of (11). (a) has the spatial symme- 
tries of pattern (a) and no spatio-temporal symmetries (cf. 2a,d); 
(b) has the symmetry properties of pattern (b) (c is one period T 
later; cf. figure 2b,e,g) 




Figure 5. Reconstructed patterns from irreps 10 and 12: 
(a,b) 10: (a, b, c) = (1, 1, 0) (cf. figure 2c,f,h); (c,d) 12: same ampli- 
tudes and same symmetries as (a,b); (e,f) 10: (a, 6, c) = (1,0,0); 
(g,h) 10: (a,6,c) = (l,l,l). 
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The second normal form (20-22) generically has three types of period-two 
points (a, 6, c): 

^(i}- 

The middle branch has the spatio-temporal symmetries of pattern (c), with 12 el- 
ements in the spatial part of the symmetry group (Ec = {kx, KyT2,Ti)). Fig- 
ure 5(a,b) illustrates this pattern (cf. figure 2c,f,h). For comparison, the pat- 
tern that would have been obtained with modes from representation 12 is in fig- 
ure 5(c,d): the symmetry group is the same, but the appearance of the pattern 
does not match the experimental observation. The first branch has a 24 element 
spatial symmetry group {p^T\, KxPTiT2,Ti) (figure 5e,f), and the third branch has 
an 18 element group (kj,T2, K^^p^, r^) (figure 5g,h). The three patterns also have 
the spatio-temporal symmetry p^tt (since = so p^ will appear in the 

symmetry group of the time-average of each of the patterns, as discussed in [1] . 

The first branch has Floquet multipliers — 1 -|- 2/x and —1 — p, (twice); 
the second branch —1 + 2/x, —1 — p+2Q ^^ ^^'^ ^" p+2Q ^^' ^^'^ ^^"^ third branch 
-1 + 2p and -1 + pq^M (twice). As a result, ii P + Q > and P + 3Q > 0, 
then all branches bifurcate supercritically, and either the first branch will be stable 
(when P < 0) or the last will be stable (when P > 0). If any branch bifurcates 
subcritically, none are stable. The middle branch, which is the one corresponding 
to the experimentally observed pattern (c), is always unstable at onset. 

6. Discussion 

Using the symmetry-based approach of Tse et al. [1] , we have analysed three exper- 
imentally observed spatial period-multiplying transitions from an initial hexagonal 
pattern. The three patterns illustrate three situations that can arise in this kind 
of analysis. Pattern (a) was straight-forward, in that a single representation of F 
had a one-dimensional space fixed by the spatial symmetry group of the pattern. 
The existence of a solution branch of the form of pattern (b) could also be inferred 
using the Equivariant Branching Lemma, though in this case it was necessary 
to include the temporal symmetry associated with period-doubling bifurcation. 
Specifically, the spatial symmetries selected a two-dimensional fixed point space 
which was further reduced to a one-dimensional fixed point space when spatio- 
temporal symmetries were taken into account. Experimentally, these two patterns 
were found for the same fiuid parameters and same 2w : 3w forcing function but 
for different frequencies w. uj = 25 Hz for (a) and a; = 35 Hz for (b). This suggests 
that the transition between these patterns, which arise for instabilities associated 
with the same representation, might be observed by tuning the frequency uj. 

Pattern (c), on the other hand, had a spatial symmetry group that fixed 
one-dimensional subspaces in two different representations, and we appealed to 
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the measured power spectrum of the pattern to choose between the two possibil- 
ities. In this situation, symmetry considerations alone were not enough. Similar 
situations arise in other bifurcation problems, for example, knowing that a stable 
axisymmetric pattern is found in a spherically symmetric bifurcation problem does 
not provide enough information to determine which is the relevant representation. 

The experimentally observed transition between hexagons and pattern (c) 
occurs by means of a propagating front that separates domains of hexagons and 
the secondary pattern. The front is initiated at the lateral boundaries of the system 
and emanates radially inward. There is little if any hysteresis, and the reverse 
transition also occurs via the same scenario. The occurrence of a front in this 
transition suggests bistability of the hexagonal pattern and pattern (c). This is 
certainly consistent with the theoretical prediction that pattern (c) is unstable at 
small amplitude, that is, at onset. However, we have not explored the possible 
stabilization mechanisms for pattern (c). 

It is worth emphasizing that an understanding of group representation the- 
ory is useful in classifying and analysing secondary instabilities of patterns, not 
only in the Faraday wave experiment as described here, but also in convection and 
other pattern formation problems (see [10]). It is also worth mentioning that the 
examples studied here indicate that spatio-temporal symmetries readily arise in 
secondary siibharmonic instabilities, and that careful experimental characteriza- 
tion of these, either by still images taken one forcing period apart or by time- 
averaging over two forcing periods, can be helpful. Subsequent instabilities of 
patterns that have spatio-temporal symmetries can be analysed using methods 
described in [11, 12]. 

The approach outlined in [1] and here is useful for taking an experimental 
observation of a secondary transition and casting it into its equivariant bifurca- 
tion theory context, but it does not predict which transitions should be expected 
in an experiment. However, in these two-frequency Faraday wave experiments, 
three- wave interactions of the type described in [13] may select a third wavevector 
that could appear in the secondary transition. Each of the representations in the 
problem under consideration is associated with a set of wavevectors, providing a 
possible mechanism for selecting between possibilities. 
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